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PREFACE 
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and  Engineering  Laboratory.  The  work  was  funded  under  DA  Project  4A762784AT42,  Design, 
Construction  and  Operations  Technology  for  Cold  Regions,  Task  BS,  Work  Unit  013,  Improve  Design 
Criteria  for  Foundations  in  Cold  Regions, 

The  author  thanks  Austin  Kovacs  and  Dr.  Devinder  Sodhi  for  technical  review  of  this  report  and 
valuable  comments. 
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Creep  and  Yield  Model  of  Ice  under  Combined  Stress 


ANATOLY  M.nSH 


INTRODUCTION 

The  successful  solution  of  ice  engineering  problems  depends  greatly  upon  the  accuracy  of  the  constitutive 
laws  and  failure  criteria  of  ice  used  in  analyses  of  engineering  structures.  Two  approaches  can  be  distinguished 
in  describing  the  deformation  and  failure  processes  of  ice  in  uniaxial  as  well  as  multiaxial  stress  states. 

In  the  first  (traditional)  approach,  it  is  assumed  that  the  total  creep  shear  strain  y®  can  be  broken  down  into 
four  components: 


f  =  Ye  +  Yp  +  Yv  +  Yt  (1) 

where  Ve,  Yp,  Yv  and  Yt  are  elastic,  primary,  secondary,  and  tertiary  pure  shear  strain,  respectively.  In  practice, 
however,  to  simplify  the  parameter  evaluation  procedure  and/or  to  obtain  closed-form  solutions  of  boundary 
problems,  some  of  the  strain  components  in  eq  1  are  often  ignored,  and  the  time-dependent  creep  process  in  ice 
is  portrayed  by  simplified  models:  primary  (yv  =  Yi = 0),  secondary  (Yp  =  Yt = 0).  and  tertiary  (Yp = Yv  =  0)  creep 
equations,  often  represented  by  various  mechanical  models.  Combined  models  such  as  attenuating  (yi  =  0), 
exponential  (Yp = 0),  and  other  creep  models  are  based  explicitly  or  implicitly  upon  various  combinations  of  the 
components  of  eq  1 .  The  failure  process,  particularly  time  to  failure,  is  either  not  considered  (in  most  cases)  or 
is  assumed  to  take  place  in  the  later,  secondary  or  tertiary  stages  of  creep.  Moreover,  often  the  ice  strength  is 
considered  a  special  field  of  study  loosely  related  to  the  creep  of  ice. 

In  the  second  approach,  deformation  and  failure  of  ice  are  considered  to  be  a  unified  process  that  takes  place 
in  all  stages  of  creep.  Consequently,  the  creep  process  is  not  considered  in  terms  of  real  time  (/)  but  in  terms 
of  nomifilized  time  t  =  where  is  time  to  failure  (Fish  1980,  Morland  and  Spring  1981,  and  others). 

Time  to  failure  plays  an  extremely  important  role  within  the  framework  of  this  approach.  It  is  one  of  the  most 
important  parameters  that  define  the  service  lifetime  of  astructure.  Time  to  failure  unites  all  three  stages  of  creep 
and  failure,  but  it  also  combines  creep  and  fracture,  microcrack  formation  kinetics,  reciystallization,  and  other 
processes  that  take  place  during  deformation  and  failure  of  ice  (Fish  1976). 

The  principal  elements  of  the  model,  which  was  originlly  developed  (Fish  1976,  1980,  1983,  1987)  to 
describe  one-dimensional  creep  and  failure  of  frozen  soils  and  ice,  and  is  expanded  below  for  three-dimensional 
defonnation,  can  also  be  approximately  expressed  in  terms  of  the  components  of  eq  1  as  a  product  of  primary, 
secondary,  and  tertiary  creep  strains,  i.e., 

f  =  Ye  +  YpYmYt  (2) 

where  Ym  =  Const.  =  viscous  shear  strain  at  failure.  Both  primary  and  tertiary  creep  strains  are  functions  of 
normalized  time,  i.e. ,  Yp  =  Yp(^  and  y,  =  ylt).  Secondary  creep  is  considered  to  be  a  point  (M)  on  the  creep  curve 
(Fig.  1 )  defining  time  to  failure.  The  stress  dependency  of  the  strain  rare  that  reaches  a  minimum  (y  =  O)  at  this 
point  is  described  by  a  viscous  flow  equation  which,  as  well  as  the  time  to  failure  function,  includes  a  yield 
criterion  of  ice. 

Thus  the  combined  creep  and  yield  model  of  ice  under  multiaxial  stress  developed  in  the  following  sections 
consists  of  four  principal  elements:  a  constitutive  equation,  a  viscous  flow  equation  and  a  yield  criterion,  all 
united  by  a  time  to  failure  function.  One  of  the  specific  features  of  this  model  is  that  any  of  its  elements  can  be 
replaced  if  it  is  found  that  other  functions  better  represent  the  mechanical  behavior  of  ice  during  creep. 


Figure  1.  Creep  model  for  ice  in  a  multiaxial  stress  state  (Fish  1983, 1984). 

CREEP  MODEL 


Triaxial  creep  tests  (ti  =  Const) 

Constitutive  equation 

The  entire  deformation  process,  which  includes  primary,  secondary,  and  tertiary  short-term  creep  of  homo¬ 
geneous  and  isotropic  ice  under  multiaxial  stress  at  constant  temperature,  can  be  presented  as  a  product  of  the 
flow  equation  and  the  nondimensional  time  function  (Fish  1980, 1983): 

(3) 

where  Yj(f)  =  =  shear  strain  rate  (intensity) 

Yjp,  =  minimum  shear  strain  rate 

he  =  -  [(si  -  £3)^  +  (^2  “  ^3)^  ~  ^1)']  ~  second  invariant  of  the  deviatoric  strain  rate 
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tensor. 

The  nondimensional  time  function 

=  1^1  =  exp{  5(r  -  In  r  -  1 )}  (4) 

where  5  is  the  first  shape  parameter  of  a  creep  curve.  Function  f(/)  is  such  that  for  primary  and  tertiar>  creep 
>  1 ,  and  for  secondary  creep  (points  M  in  Fig.  1 )  r = /m>  t  =  1 .  and  1 .  For  ice,  the  value  of  5  is  between 

0.3  and  0.7  and  has  an  average  value  of  -0.5  (Fish  1987). 
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Flow  equation 

The  minimum  shear  strain  rate  is  related  to  the  time  to  failure  by  means  of  the  flow  equation 

Yini  =  —  or  Yin,fm  =  C  (5) 

where  C  =  the  second  shape  parameter.^hich  is  equal  to  the  viscous  strain  at  failure.  For  polyciystalline  ice 
C  is  between  10"'  and  10"^Both5andC  are  assumed  to  be  constant  for  a  given  ice  temperature.  The  tilde  (~) 
indicates  that  in  the  general  case  C  may  be  a  function  of  stress,  temperature  and  other  factors  (Fish  1987, 1 989). 


Time  to  failure 

When  time  t  =  the  strain  rate  reaches  a  minimumYi  =  iim  (y  j  =  O),  and  it  is  assumed  that  failure  occurs. 
In  the  framework  of  this  model,  time  to  failure  tm  performs  a  crucial  role:  it  not  only  unites  all  stages  of  creep, 
which  makes  it  possible  to  describe  the  entire  creep  process  under  both  constant  stresses  and  constant  strain  rates 
by  means  of  a  single  constitutive  equation,  but  it  also  combines  the  long-term  strength  and  the  yield  criteria  of 
ice  so  that  the  latter  can  be  incorporated  into  the  constitutive  equation  of  ice.  The  magnitude  of  shear  stress  that 
causes  failure  at  time  t^  can  be  expressed  (Fish  1991)  as  a  product  of  two  independent  functions:  the  yield 
function  (criterion)  and  the  nondimensional  time  function  <i)(^: 

(6) 

where  J20  =  J2o(^i )  is  an  “instantaneous”  yield  criterion,  I\  is  the  first  invariant  of  the  stress  tensor,  and.f2  is  the 
second  invariant  of  deviatoric  stress  tensor.  By  selecting  the  time  function  in  the  simplest  form, 

where  t  =  /m/to  =  nondimensional  time  to  failure 
/o  =  relaxation  time  (see  eq  12  below) 
n  =  parameter. 


Combining  eq  6  and  7,  the  failure  or  the  long-term  strength  criterion  of  ice  takes  the  form 


(8) 

(9) 


where  Tj  =  and  =  7 can  be  defined  as  the  shear  stress  and  the  “instantaneous”  shear  stress  intensity 
(resultant)  respectively  (see  eq  38  below). 


Modified  flow  equation 

Combining  eq  5, 8  and  9,  the  flow  equation  of  ice  in  a  multiaxial  stress  state  becomes 

i«/2 


V:.  = 


.c  iL 

/q  '’^io 


(10) 


or 


L) 


(11) 
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where 


and  to  =  Frenkel’s  relaxation  time  (Frenkel  1947) 
h  =  Planck’s  constant 
k  =  Boltzmann’s  constant 
T  =  absolute  temperature 
E  =  activation  energy 
/?  =  gas  constant. 

Note  that  the  flow  equations  (eq  10  and  1 1)  are  fundamentally  different  from  th^Norton-Glen  (1958)  flow  law 
(see  eq  43  below)  although  it  contains  a  power  function  of  stress.  Parameters  C  and  to  in  eq  10  wd  1 1  have  a 
deflnite  physical  meaning,  and  the  denominator  of  the  stress  function  is  a  temperature'dependent  yield  criterion 
of  ice  (see  eq  38  below).  Hiis  yield  criterion,  which  is  a  function  of  the  first  invari  wt  of  the  stress  tensor,  relate:' 
the  minimum  shear  strain  rate  and  the  shear  stress  in  the  whole  spectrum  of  hydrostatic  pressures  (mean  norm  i 
stresses).  At  the  mean  stress  Om  <  Omax  the  first  invariant  serves  as  a  depressant  forthe  strain  rate;  i.e.  the  increa  se 
of  the  hydrostatic  '.tress  is  accompanied  by  a  decrease  of  the  strain  rate.  When  stress  On,  >  Omax  the  increase  of 
the  hydrostatic  str  iss  enhances  the  creep  rate  that  is  confirmed  by  test  data  (Jones  and  Chen  1983), 

The  original  Norton  equation  is  more  suitable  for  frictionless  materials,  such  as  metals,  the  creep  behavior  of 
which  is  unaffected  by  the  first  invariant  of  the  stress  tensor.  If  the  combination  of  applied  stresses  is  such  that 
Ti  =  Tjo  and  rm  =  /q  ihen  eq  10  becomes 

Yio=^  (13) 

^0 

where  Yj^  =  Yj^  can  be  defined  as  the  “instantaneous”  shear  strain  rate,  i.e.,  the  rate  at  which  the  transition  from 
the  brittle  to  the  ductile  mode  of  failure  takes  place. 


Creep  shear  strain 

At  constant  stresses  creep  shear  strains  are  obtained  by  integrating  eq  3.  The  entire  creep  process  as  described 
by  (Fish  1983)  is 


T^(f)  =  l+Yir\|/ 
G 


where  Tj  =  shear  stress 
G  =  shear  modulus 
Yf  =  2J20  =  creep  shear  strain 


•^26  =  -  [(ci  -  hf  +  “  ^3)^  +  (^3  ~  ^1)^]  =  second  invariant  of  the  deviatoric  strain  tensor 
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V  =  =  integration  coefficient  calculated  depending  upon  the  normalized  time  t . 


(14) 


In  the  practical  sense,  the  instantaneous  strains  are  small  and  in  many  instances  can  be  ignored,  which 
substantially  simplifies  the  parameter  evaluation  procedure.  When  t  =  tm,t=\,  and  yj  =yj^,  eq  14  gives  an 
expression  for  failure  shear  strain  (Fig.  1) 


Yfm  J  Yim^mVm  (1^^ 

where  Vm  is  an  integration  coefficient  calculated  for  r  =  1.  The  first  term  in  the  right  side  of  eq  15  is  the 
instantaneous  strain  and  is  ignored  in  the  plot  shown  in  Figure  1. 

Combining  eq  5  and  15 


(16) 
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where  «  1/(1  -  8)’^  =  Md  X  =1  — 5..it  is  more  cohyenjeht  fofthe^nalysis  bftest  data-to.u^^^^  . 
approximate  integral  of  eq  3  (Fish  1983,  eq  38;  Fish  1987,  eq  20): 

Y^(<)  “  ^  +  -pr  ?  V‘^'~  (17) 

n 

where  Sj  =  X^^-X.  Defining  C(r)  =  and  combining  eq  5  and  17  the  entire  creep  process  can  be  expressed 

in  terms  of  viscous  (nonrecoverable)  shear  strain  only: 

if (r)  =  3.  +  -£ (18) 

^  ix\cj 


where  82  =  5  i/C  . 

Thus,  the  entire  cr^  process  in  ice  in  a  multiaxial  stress  state  at  constant  stresses  is  described  by  means  of 
only  two  parameters,  C  and  S{X),  to  be  determined  from  test  data.  An  evaluation  procedure  for  these  parameters 
is  presented  in  Fish  (1983, 1987, 1989). 


Triaxial  constant  strain  rate  tests  (yj  =  Const.) 


Strength  criterion 

If  the  tests  are  carried  out  in  such  a  manner  that  the  shear  strain  rate  is  maintained  constant  (strength  tests) 
at  Yj  =  Yjn,,  the  shear  stress  becomes  a  function  of  this  rate.  Defining  Yj  =  Yjt.  Yim  “  Yim  W  ond  C  =  Yin,i  and 
combining  the  latter  with  eq  10, 


Mm 


Mo 


(19) 


where  tjo  is  a  yield  criterion  (eq  38  below,).  Subscript  m  indicates  that  the  magnitudes  of  stress  (Tjm)  and  strain 
(Yim)  are  referred  to  the  point  M  corresponding  to  the  maximum  (peak)  shear  stress  in  Figure  1. 

If  Yj  =  Yio  Yim  “  Yio  ®  Yio^o  ''im  "  "'io*  a™l  eq  19  will  coincide  with  eq  1 3;  i.e.,  the  maximum  (peak) 
shear  stress  will  coincide  with  the  yield  stress. 


Stress-strain  relationship 


Stress-strain  relationships  outside  point  M  (Fig.  1)  are  obtained  by  combining  and  rearranging  eq  3, 4,  and 


where  t  =t/t^  =  Yi/Yjn,  =  Y.  and  Tj  =  Ti(^  =  Ti(Y).  For  secondary  creep,  r  =  1  and  the  strain  rate  dependency  of 
stress  is  defined  by  eq  19. 


YIELD  MODEL 
Extended  von  Mises  criterion 

Analysis  of  test  data  shows  that  at  low  stress  levels  (Haynes  1975)  or  low  strain  rates  (Hawkes  and  Mellor 
1972)  the  strength  of  ice  is  independent  of  confining  pressure,  the  angle  of  internal  friction  of  ice  approaches 
zero,  and  ice  may  be  considered  as  an  ideally  cohesive  material  whose  yield  criterion  may  be  selected  in  the  form 
of  either  the  von  Mises  or  the  Tresca  rupture  models.  On  the  other  hand,  triaxial  tests  carried  out  at  relatively 
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high  stresses,  high  strain  rates  and/or  high  confining  pressures  (Sayles  1974,  Jones  1978,1982,  Hausler  1982, 
Richter-Menge  et  al.  1986  and  others)  indicate  that  the  strength  of  ice  is  a  nonlinear  function  of  the  confining^ 
pressure.  Sayles  ( 1 974)  showed  that  the  angle  of  internal  friction  of  ice  differs  from  zero.  Moreover,  at  a  certain 
level  of  confining  pressure  the  strength  of  ice  reaches  a  maximum  value  (Jones  1978)  associated  with  the  ice 
melting  pressure.  The  yield  function  of  such  a  material  can  be  expressed  in  terms  of  the  extended  von  Mises 
or  Mohr-Coulomb  rupture  models  (Fish  1991),  and  in  the  former  case  by  expanding  the  vori  Mises  criterion  into 
the  series 

j]F  =  d+  </,/,  +  +  d^l]  +  d/^  . . .  (22) 

where  J2  =  second  invariant  of  the  deviatoric  stress  tensor 
1 1  =  first  invariant  of  the  stress  tensor 
^.^1, 2,3,4  =  parameters. 

By  retaining  only  the  first  term  on  the  right  side  of  eq  22  it  becomes  von  Mises’  yield  criterion:  the  first  two 
terms — the  von  Mises-Drucker-Pragercriterion;  three  terms — Nadai’s  type  yield  criterion  (Nadai  1 950;  Smith 
1 974;  and  others);  four  terms — the  teardrop  model  for  ice  (Nadreau  and  Michel  1 986, 1 986a);  etc.  We  will  limit 
ourselves  to  consideration  of  the  second-order  polynomial  function  as  the  simplest  nonlinear  yield  function 
whose  parameters  can  be  easily  determined  from  test  data. 

Parabolic  yield  criterion 

Let  us  assume  that  the  compressive  stresses  are  positive.  Then  the  strength  of  ice  in  a  multjaxial  stress  state 
can  be  described  by  a  parabolic  yield  criterion  in  Figure  2  (Fish  1991): 

Ti  =  c+ha„,-:^o^  (23) 

where 

tj  =  4^-  =  ^[(o,-02)2  +  (02-03)2  +  (03-0,)2]'^ 
is  shear  stress  (intensity)  and 

On,  =  —  =  -  (O]  't-  02  +  03)  =  mean  normal  stress 


2 


Figure  2.  Yield  criteria  of  ice:  I -von  Mises,  2-Drucker-Prager,  3-Paraholic  (Fish  1991 ). 
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oi?  02, 03=  principal  stresses 

c  =  resistance  in  pure  shear  associated  with  the  ice  cohesion  on  the  octahedralpl^e 
b  =  tanij)  and  (])  =  angle  of.intemal  friction  on  the  octah^ral  plane 


Omax  =  ics  melting  pressure,  i.e.,  the  magnitude  of  the  mean  normal  stress  at  which  the  shear  stress  of  ice 
reaches  a  maximum: 


maxTi  =  c+|o^a^. 

The  yield  curve  intersects  the  hydrostatic  axis  at  points  Hi, 2,  the  abscissas  of  which  are 
IH,2=  ‘^max  ±  [<^ax  +  ^  Omax]’^- 


(24) 

(25) 


Note  that  the  yield  curve  shape  is  defined  by  the  ratio  of  parameters  c  and  b. 

Equation  23  can  be  considered  as  an  extended  von  Mises-Drucker-Prager  yield  criterion.  At  low  stress  levels 
(®max  °°)>  transforms  into  the  Drucker-Prager  (1952)  yield  criterion.  For  frictionless  materials  (b = 0), 

eq  23  reduces  to  the  von  Mises  yield  criterion.  Parameter  Cmax  niay  be  regarded  as  one  of  the  fundamental 
physical  characteristics  of  ice,  probably  closely  related  to  thfe  phase  changing  pressure  (point  Hi)  for  a  given 
temperature  (Nadreau  and  Michel  1986;  Hallam  and  Nadreau  1988).  Point  H2  with  abscissa  h2  may  be 
considered  as  the  resistance  of  ice  to  triaxial  extension.  Introduction  of  a  new  ice  physical  characteristic  Opiax 
and  ratio  fc/2amax  23  makes  the  latter  different  from  the  known  parabolic  yield  criteria  of  Nadai  (1950), 

Smith  (1974),  and  others.  Instead  of  the  ratio  i>/2CTniax  these  authors  use  the  parameter  ^2  (eq  22),  the  physical 
meaning  of  which  is  uncertain  and  the  magnitude  of  which  will  be  dependent  of  the  strain  rate  (or  time  to  failure). 
Therefore,  these  criteria  cannot  be  applied  to  describe  the  test  data  plotted  below  with  a  single  set  of  the  strength 
parameters. 

Thus  ice  strength  in  a  multiaxial  stress  state  can  be  described  by  a  yield  criterion  (eq  23),  tnree  parameters 
of  which — c,  b(^)  and  Omax — have  a  definite  physical  meaning  and  can  be  determined  from  test  flata. 


LONG-TERM  STRENGTH 


Triaxial  creep  tests  (q  =  Const.) 

It  has  been  mentioned  above  that  the  creep  strength  of  ice  is  a  function  of  time  or  the  applied  strain  rate.  The 
absolute  values  of  ice  strength  as  well  as  the  magnitudes  of  the  strength  parameters  change  from  the  greatest 
(instantaneous)  values  at  the  time  of  loading  (fg)  to  zero  for  longer  time  intervals.  Consequently  the  creep 
strength  of  ice  will  not  be  characterized  by  a  single  yield  curve  but  by  a  family  of  curves  for  given  times  (g  < 
(Fig.  3)  or  given  strain  rates  Yig  >  Yj  j  >  Yj^ .  During  creep  at  constant  stress,  the  strength  parameters  in  eq 
23,  except  Oniax,  become  functions  of  time  to  failure: 

Ti(0  =  c(f)  +  b(f)  o2  (26) 

^^max 

where  c(f)  and  b(/)  =  tan$(f)  are  the  cohesion  and  this  friction  angle,  on  the  octahedral  plane,  at  time  /  =  fm, 
respectively.  The  ice  strength  reaches  its  maximum  w  hen 

maxTi(r)  =  c(r)+^On,ax-  (27) 

Note  that  the  shape  of  the  creep  strength  curves  and  the  abscissas  of  points  Hi  and  H2  in  Figure  3  are  dependent 
upon  the  ratio  of  parameters  c(/)  and  b(t)  and  their  change  with  time: 


(28) 
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Ti  =  V^ 


0^^  =  -T-  =  Mean  Normal  Stress 
o 


t,(t)  =  c(t)  +  b(t)dni- 


b(t). 

20|riax 


max  Tj  (t)  =  c  (t)  + 


2 


Omax 


bl.2  =  Cmax* 


2c  (t) 
b(l) 


b(t)  =  tan(|i(t) 


cW 

b{t) 


=  Const. 


Figure  3.  Creep  strength  criterion  (Fish  199J), 


As  a  first  approximation  one  may  assume  that  this  ratio 

^  =  Const  (29) 

bif) 

i.e„  the  coordinates  of  the  points  H],2  during  creep,  remain  unchanged  and  all  the  curves  will  have  common 
intercepts.  At  low  stress  levels  or  when  o^ax  -♦  eq  26  reduces  to  the  von  Mises-Drucker-Prager  strength 
criterion  expanded  for  creep  conditions: 


Ti(r)  =  c(r)  +  6(r)on,  (30) 

and  when  6  =  0  it  reduces  to  the  von  Mises  strength  criterion 

Ti(0=c(f),  (31) 

in  which  parameters  c  and  b  are  functions  of  time. 

Triaxial  constant  strain  rate  tests  (y,  =  Const  .) 

For  constant  strain  rate  (strength)  tests  =  Y/Yi  =  Yin/Yim  >  “  Yio/Yjo  •  parameters  of  eq  26 

through  31  as  well  as  the  strength  criteria  of  ice  become  functions  of  the  shear  strain  rate  (and  strain),  i.e.: 


ti(Yi)  =  c(Yi)  +  6(Yi)om  “ 

^^max 

(32) 

maxTi(Yi)-c(Yi)+ 

(33) 

r  5  2c(y)  1*^ 

^l,2(Yi) ~  ®max  ^  ^ax  ®max 

(34) 

When  Omax  oo  and  b*0, 

'ti(Yi)=c(Yi)+ 6{Yi)omax-  (35) 
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When  6  =  0, 


(36) 


For  constant  strain  rate  tests  the  creep  strength  criterion  depicted  in  Figure  3  will  retain  its  original  meaning. 
However,  the  time  functions  of  the  strength  parameters  in  it  must  be  replaced  by  functions  of  shear  strain  rates 
(Fig.  4).  The  creep  strength  curves  will  have  common  intercepts — ^points  if  the  ratio 


=  Const. 


(37) 


If  this  ratio  varies,  i.e.,  if  the  rates  of  change  of  parameters  c  and  b  are  different,  the  curves  will  intersect. 


Time-dependent  failure 

Passing  vertical  planes  through  point  M  in  Figure  3  one  obtains  the  relationship  between  the  creep  strength 
of  ice  Tj  and  the  time  to  failure  /m.  Such  a  dependency  relates  the  (instantaneous)  yield  function  of  ice  determined 
from  creep  tests  at  time  to  (or  from  constant  strain  rate  tests  at  the  shear  strain  rate  Yjo)  with  a  err  ep  strength  Tjf/) 
=  determined  at  times  /m  =  ^l<  clc- 

It  has  been  shown  elsewhere  (Fish  1991)  that  the  creep  strength  criterion  (eq  26)  can  be  expressed  as  a  product 
of  the  yield  function  (eq  23)  and  the  nondimensional  time  function  (eq  7),  i.e. 


(6) 

when  /m  =  to  >  ~  ^  ^nd 

Tjo  =  c  +  6oOp,  — 

(38) 

^0  “  •  ^0  “  tan<|)Q. 


Subscript  0  indicates  that  the  failure  shear  stress  or  the  strength  parameters  are  referred  to  the  instantaneous 
yield  conditions.  Comparing  eq  6, 7, 9, 26,  and  38 


Consequently,  for  constant  strain  rate  tests  (tj  =  Xjm)  eq  39  takes  the  form 

=  =  =  (40) 

'10  K 

Note  that  the  two  parameters  to  and  n  define  the  shapes  of  the  curves  of  the  long-term  strength  of  ice.  Analysis 
of  Jones’  data  (Nadreau  and  Michel  1986a)  showed  that  for  ice  under  triaxial  compression  the  exponent  n 
changes  insignificantly.  However,  the  magnitude  of  n  increases  slightly  when  the  confining  pressure  is  zero. 
It  has  been  shown  (Fish  1 99 1 )  for  frozen  soils  at  six  different  types  of  loading  (uniaxial  compression  and  tension, 
pure  shear,  and  triaxial  compression  at  three  different  mean  normal  stresses)  that,  for  practical  purposes,  these 
parameters  may  be  considered  to  be  independent  of  the  loading  regime. 

Thus  the  entire  process  of  defonnation  and  failure  of  ice  in  a  multiaxial  stress  state  at  both  constant  stresses 
and  constant  strain  rates  is  described  by  seven  parameters: 

Creep  5(k)andC 
Failure  time  /oandn 

Yield  Co,  60(90).  and  CTmax 
which  are  determined  from  test  data. 
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On,,  Mean  Normal  Stress  (MPa) 


Figure  4.  Strength  of  polycrystalline  ice  under  triaxial  compression  (aj>  02 
=  03)  at  -11.8°C  (data  from  Jones,  1982). 

TEST  DATA 

Yield  of  ice  under  triaxial  compression 

The  triaxial  constant  strain  rate  tests  used  in  the  present  analysis  were  performed  by  Jones  (1982).  The 
freshwater  ice  used  in  them  had  a  randomly  oriented  polycrystalline  structure.Tlie  laboratory-made  ice  samples 
were  20  mm  in  diameter  and  60  mm  long.  The  axial  strain  rate  e|  varied  between  10"^  and  10“’  s“’  and  the 
confining  pressure  (02  =  03)  ranged  from  0. 1  to  85  MPa.  The  tests  were  carried  out  at  -1 1 .8°  ±0.9®C.  The  test 
data  were  recalculated  by  this  author  in  terms  of  the  shear  stress  (Xj)  versus  the  mean  normal  stress,  Om,  and  are 
presented  in  Figure  4. 

Strictly  speaking,  this  type  of  triaxial  test  can  only  approximately  be  considered  as  a  pure  constant  strain  rate 
regime.  Indeed,  during  each  test  the  axial  strain  was  applied  at  a  certain  rate  (constant  strain  rate  regime),  while 
the  confining  pressure  which  corresponds  to  a  creep  regime  at  constant  stress  was  maintained  constant. 

Nevertheless,  as  a  first  approximation  one  may  assume  that  eq  19  can  be  applied  to  this  type  of  test.  For  tests 
when  O]  >  02  =  03  it  takes  the  form 


Jones  also  found  that  exponent  n  =  3.95,  and  that  ai  strain  rates  ej  =  5.4  x  10“?  s“*  and  above  the  ice  strength 
practically  does  not  change.  Therefore,  for  this  particular  type  of  ice  and  for  the  given  test  temperature  this  rate 
can  be  considered  an  “instantaneous”  strain  rate.  Since  eq  42  and  43  have  the  same  (power)  structure  the 
exponent  n  can  be  used  in  the  following  analysis. 

Then  the  magnitudes  of  the  instantaneous  shear  stress  as  a  function  of  the  mean  stress  tjo  =Aom)  can  easily 
be  calculated  by  eq  42  using  test  data  from  Figure  4.  After  that  the  parameters  of  the  yield  criterion  can  be 
calculated  by  eq  38  and  25. 

Based  upon  a  regression  analysis  of  all  86  tests  the  parameters  of  the  yield  criterion  turned  out  to  be 
Co  =  16.4  MPa 
bo  =  0.27 
<1)0  “  15° 

c^max  “41.5  MPa 
hi  =  123.7  MPa 
/i2  = -40.7  MPa. 

Both  the  experimental  and  predicted  (by  eq  38,  40  and  42)  values  of  the  ice  shear  strength  under  triaxial 
compression  at  various  axial  strain  rates  and  at  various  mean  normal  stresses  are  presented  in  Figure  4.  One  can 
conclude  that  equations  developed  fit  the  test  data  quite  well. 

It  is  interesting  to  note  that  the  magnitude  of  Omax  (41-5  MPa)  is  in  agreement  with  that  found  by  Jones 
(-40  MPa).  The  magnitude  of  the  intersect  (hi  =  123.7  MPa)  is  in  agreement  with  the  value  (-120  MPa) 
suggested  in  Hallam  and  Nadreau  (1988). 

Note  that  the  instantaneous  value  of  the  ice  friction  angle  turned  out  to  be  very  small  ((l)o  =  15°).  Both  the 
friction  angle  and  the  cohesion  decrease  rapidly  with  decreasing  strain  rate,  and  at  a  strain  rate  of-lO"^  s"'  the 
friction  angle  magnitude  <1) « 1  °.  That  means  that  the  ice  gradually  transforms  from  a  material  that  at  high  strain 
rates  possesses  both  the  cohesion  and  the  friction  angle  into  an  ideally  cohesive  material  the  shear  strength  of 
which  may  be  described  by  a  the  von  Mises  (eq  36)  or  the  Tresca  strength  criterian  expanded  for  creep  con¬ 
ditions.  The  latter  is  a  particular  case  of  a  more  general  parabolic  (extended)  Mohr-Coulomb  strength  criterion 
that  also  can  be  applied  to  describe  the  entire  family  of  yield  surfaces  in  Figure  4  (Fish  1991). 

Creep  of  ice  under  triaxial  compression 

Triaxial  tests  of  ice  used  in  the  creep  analysis  were  performed  by  Golubov  et  al.  (1990)  using  laboratory- 
made  samples  of  saline  ice.  The  test  samples  were  45. 1  mm  in  diameter  and  100  mm  long.  The  diameter  of  the 
randomly  oriented  ice  crystals  varied  between  0.5  and  1 .5  mm.  The  density  of  the  ice  was  0.85  to  0.9  g/cm^  and 
average  salinity  was  -5%o.  The  tests  were  carried  out  at  -5°C  ±0.2°C. 

Prior  to  the  creep  tests  a  series  of  strength  tests  were  conducted  at  a  constant  strain  rate  and  various  confining 
pressures  (02 =03).  The  data  from  these  tests  were  recalculated  in  terms  of  tj  and  are  presented  in  Figure 

5 .  The  strain  rate  at  which  the  tests  were  conducted  was  not  specified,  but  comparison  with  the  creep  tests  shown 
in  Figure  6  brings  us  to  the  conclusion  that  the  strength  test  strain  rate  was  higher  than  the  secondary  creep  rate 
in  the  creep  tests.  Therefore,  the  data  in  Figure  5  can  be  used  for  determining  the  strength  parameters  of  the  saline 
ice.  The  tests  were  conducted  under  considerably  smaller  mean  normal  stresses  than  the  tests  of  Jones  (Fig.  4). 
Therefore  the  test  data  can  be  approximated  by  the  linear  von  Mises-  Drucker-Prager  yield  criterion  (eq  38  at 
fifo  =  0).  The  magnitudes  of  the  strength  parameters  of  the  saline  ice  at  -  5°C  turned  out  to  be 

Co  =  0.58  MPa  and  &o  =  0.11  (())o  =  6°) 

which  are  considerably  smaller  than  those  of  freshwater  polycrystalline  ice  at  -1 1 .8°C. 

The  results  of  the  triaxial  creep  tests  are  presented  in  Figure  6.  These  tests  were  carried  out  at  constant  shear 
stress  Tj  =  (o]  -  03)  /  VJ  =  0.46  MPa  and  at  three  different  magnitudes  of  the  mean  normal  stress  (am  =1.18, 
2. 1 8  and  4. 1 8  MPa).  From  a  cursory  analysis  of  the  test  data  ( ne  can  conclude  that  when  all  other  conditions 
are  equal  the  creep  shear  strain  and  the  strain  rate  in  ice  substantially  depend  upon  the  mean  normal  stress  as 
predicted  by  eq  10. 
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Figure  5.  Strength  of  saline  ice  under  triaxial  compression  (o/  >  02  =  cj)  at 
-5°C  (data  from  Golubov  et  al,  1990), 


Figure  6.  Creep  of  saline  ice  at -5  °C  under  triaxial  compression  (a  i  >  02 = aj)  at  constant 
shear  stress  x,-  =  0.46  MPa  and  various  mean  normal  stresses  <5^  =  1 .18;  2.18  and  4.18 
MPa  (data  from  Golubov  et  al.  1990). 
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The  plot  in  Figure  6  contains  both  the  experimental  and  predicted  creep  curve^  The  latter  were  calculated 
by  eq  5, 9, 1 0, 1 7  and  38  using  the  following  magnitudes  of  the  creep  parameters:  C  =  0. 1 ;  X = 0.57;  6i  =  0. 1 85; 
n  =  5.46;  to  =  6.63  x  10^  s;  and  cq  and  bo  obtained  from  the  test  data,  assuming  that  the  instantaneous  strains 
are  small  and  can  be  ignored.  One  can  conclude  that  correspondence  between  the  experimental  and  predicted 
creep  shear  strains  is  quite  good,  taking  into  account  the  fact  that  a  single  set  of  parameters  was  used  in  the 
calculations,  together  with  a  very  limited  amount  of  test  data. 


CONCLUSIONS 

1 .  A  combined  creep  and  yield  model  has  been  developed  for  ice  under  multiaxial  stresses  for  both  constant 
stress  and  constant  strain  rate  loading  regimes.  The  creep  model  (Fish  1 980, 1 983)  describes  the  entire  (primary, 
secondary  and  tertiary)  creep  process  in  ice  in  terms  of  shear  strain  and  consists  of  four  principal  elements:  a 
constitutive  equation,  a  flow  (secondary  creep)  equation  and  a  yield  criterion,  all  united  by  a  time  to  failure 
function. 

2.  Secondary  creep  is  considered  to  be  a  transition  point  from  primary  to  a  tertiary  creep  stages  defining  time 
to  failure.  The  strain  rate  at  this  point  is  a  function  of  both  the  second  invariant  of  the  deviatoric  stress  tensor 
(^2)  (shear  stress)  and  the  first  invariant  of  the  stress  tensor  (/  j )  (hydrostatic  pressure).  The  stress  dependency 
of  the  secondary  (minimum)  creep  strain  rate  at  failure  is  described  by  a  modified  Norton-Glen  (power)  equa¬ 
tion,  which  as  well  as  the  time  to  failure  function,  include  a  yield  criterion  for  ice. 

3.  The  yield  model  is  selected  in  the  form  of  the  (parabolic)  extended  von  Mises-Drucker-Prager  or  the  ex¬ 
tended  Mohr-Coulomb  yield  criterion  (Fish  1 99 1 )  in  which  the  ice  strength  is  characterized  by  three  parameters: 
the  cohesion  c  and  the  friction  angle  (|»,  on  the  octahedral  plane,  and  the  ice  melting  pressure  Omax.  the 
magnitude  of  the  mean  normal  stress  at  which  the  shear  strength  of  ice  reaches  a  maximum  value. 

4.  The  model  has  been  verified  using  test  data  on  the  strength  of  freshwater  polycrystalline  ice  at  -1 1 .8°C 
(Jones  1982)  and  creep  of  .saline  ice  at  -5°C  (Golubov  et  al.  1990),  both  under  triaxial  compression. 
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